Electrostatic Field Driven Alignment of Organic Oligomers on ZnO Surfaces 
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We study the physisorption of organic oligomers on tiie ZnO (1010) surface using first-principles 
density-functional theory and non-empirical embedding methods. We find that both in-plane loca- 
tion and orientation of the molecules are completely determined by the coupling of their quadrupole 
moments to the periodic dipolar electric field present at the semiconductor surface. The adsorption 
is associated with the formation of a molecular dipole moment perpendicular to the surface, which 
bears an unexpected linear relation to the molecule-substrate interaction energy. Long oligomers 
such as sexiphenyl become well-aligned with stabilization energies of several 100 meV along rows of 
positive electric field, in full agreement with recent experiments. These findings define a new route 
towards the realization of highly-ordered self- assembled arrays of oligomers/polymers on ZnO (1010) 
and similar surfaces. 

PACS numbers: Valid PACS appear here 
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Hybrid structures made of conjugated organic 
molecules and inorganic semiconductors exhibit an enor- 
mous application potential as they combine the favorable 
features of both components in a single new material 1] . 
However, interfacing of organic molecules with the typ- 
ically highly reactive semiconductor is a complex issue. 
Rupture and fragmentation are frequently observed lead- 
ing to ill-defined interfaces [2}. On the other hand, the 
electronic structure of the semiconductor surface might 
be exploited for developing novel strategies of molecular 
aggregation. In this Letter, we demonstrate that the elec- 
trostatic interaction between the semiconductor and the 
TT-electron system gives indeed rise to the self-assemblage 
of stable and highly ordered monolayers for a wide class 
of conjugated organic molecules. 

The specific surface under consideration is the non- 
polar (1010) crystal plane of ZnO. The chemistry of ZnO 
surfaces, see e.g. Ref. i, has been largely investigated 
in the context of catalysis [3] and, more recently, much 
attention is paid to the linkage with organic dyes and 
polymers, driven, e.g., by photovoltaic applications [^Q. 
In particular, it has been found experimentally that p- 
sexiphenyl (6P) absorbs flat on the ZnO(lOTO) surface 
with the long axis of the molecule perpendicular to the 
polar [0001] direction [3]. In this study, the hybrid in- 
terface has been formed entirely under ultra-high vac- 
uum conditions suggesting that intrinsic features of the 
semiconductor-molecule system are behind that type of 
aggregation. The theoretical analysis presented below 
not only confirms this conjecture but reveals systematic 
tendencies common to all oligomers that can be used to 
engineer the growth of inorganic/organic structures. 

In order to establish a proper and efficient methodical 
framework, we start with biphenyl (2P) as a short model 



oligomer. Fig. [T^ and b depict the configuration ex- 
amined. The origin of the reference coordinate system is 
located at the center of a surface Zn-0 bond, the z- and y- 
axis point along the surface normal and the polar [0001] 
direction, respectively. The position of the molecule is 
denoted by the coordinates of its center of mass. We con- 
sider a clean, non reconstructed surface optimized using 
density functional theory (DFT) as described in Ref. S 
Our goal is the construction of the ground-state potential 
energy surface (PES) of the molecule-semiconductor sys- 
tem. In a first step, we set the center of the molecule on 
top of the Zn-0 bond {x~y^O) with its long axis aligned 
in x-direction and the molecular plane parallel to the sur- 
face. The interaction energy of this arrangement keeping 
both the molecular and ZnO(lOlO) surface configuration 
frozen is plotted versus distance from surface (z) in Fig. 
[TJ;. The curves are computed at the PBE level [9| and 
with dispersion correction (PBE-I-D) [lO|, using two dif- 
ferent computational methods: a periodic pseudopoten- 
tial plane- wave (PW) approach[lli and the periodic elec- 
trostatic embedded cluster method (PEECM) [H, 12|, as 
implemented in the TURBOMOLE ^13] program. 

Figure [1]; shows that the PEECM results agree very 
well with the PW ones. The practical advantage of 
PEECM lies in the fact that it considers only a single 
molecule in interaction with the surface. Unlike the PW 
method, where a whole periodic organic monolayer (of 
hypothetical structure) has to be treated, PEECM de- 
fines thus a cost effective way to tackle the initial ad- 
sorption step of the molecule. As expected from previ- 
ous studies, see e.g. Ref. [l4|, the PBE functional leads to 
weak binding, while the dispersion correction increases 
the binding energy and reduces the molecule-substrate 
distance. In order to verify the accuracy of PBE+D for 
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FIG. 1: a) Top and b) side view of 2P on the 
ZnO( 1010) surface (dotted rectangle: surface unit cell), c) 
Molecule-substrate interaction energy for 2P verus vertical 
distance z {x='y—0), as computed by different theoretical 
methods (see text). 



ZnO surfaces, we performed reference MP2 calculations 
within the PEECM scheme [ll|. Figure [T]; exposes that 
PBE-I-D is quite far from MP2 and thus cannot be safely 
used for ZnO surfaces. The MP2 predicts an interaction 
energy of 370 meV with an equilibrium distance zq w 3.5 
A. This molecule-substrate distance will be used in all 
the following calculations. 
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FIG. 2: a) Molecule-substrate interaction energy versus y- 
position measured relative to y = (x = 0, 2 = zo), by differ- 
ent theoretical methods. The length of the ZnO(lOlO) surface 
unit cell in this direction is 5.19 A. b) Induced molecular 
dipole moment, components /iy and (nx = 0) computed 
in QM/PPC. The electrostatic field of the ZnO (1010) surface, 
from PW/PBE calculations, is illustrated (color online) in 
c) (^-component) and d) (j/-component). The colormap cov- 
ers the range from -5V/nm (intense blue) to 5V/nm (intense 
red). Also shown is the molecule position at minimum Hy 
(2/ = 2.5 A). 



As a next step, we now consider the change of the 
molecule-substrate interaction energy (AE) when the 
molecule is translated along y-direction. As displayed in 
Fig. the computations performed again in different 
approximations commonly reveal a distinct minimum for 
y « 3.7 A, i.e., when the center of the molecule is close to 
a position atop a Zn atom. Interestingly, the relative en- 
ergies of the PBE(PW) approach are practically not mod- 
ified by inclusion of the dispersion. In order to save com- 
putation time for the examination of the complete PES 
below, we performed the same procedure but modeling 
the ZnO(lOlO) surface only by point-charges with values 
+q and —q at the lattice positions of the Zn and O atoms, 
respectively. We call this method QM/PPC (quantum 
mechanics/periodic point-charges), because the molecule 
is treated quantum mechanically at PBE level, while the 
ZnO surface is classically described. Hence, only the elec- 
trostatic interaction between the molecule and the sub- 
strate is considered in this approach. For q—1.2, excellent 
agreement with the PBE(PW) result is indeed achieved. 
This value q is very close to what is found in the MuUiken 
population analysis of the ZnO(lOTO) surface Q. There- 
fore, we conclude that exchange-correlation forces deter- 
mine the absolute energy (see Fig. Ic), but the energy 
variation when moving the molecule within the surface 
plane is completely dominated by the electrostatic cou- 
pling. Exchange-correlation effects vary on the atomic 
length scale, but are averaged out as the molecule is 
larger than the ZnO(lOlO) unit cell. 

The alternating point charges which characterize the 
ZnO(lOlO) surface create a periodic dipolar electric field 
F. An important consequence of this field is that it gen- 
erates in turn an induced dipole moment fl in the 2P 
molecule. For symmetry reasons, is negligible, while 
Fy and F^ reach values of several V/nm. The resultant 
fly and Hz are plotted versus y-position in Fig. ^jp. Over 
the length a of the unit cell, they change sign with a rel- 
ative of shift of a/4. This behavior reflects the dipolar 
character of F, as illustrated in Fig. [2]:,d. The electric 
field is largely inhomogeneous, but sufficiently far from 
the surface where the molecule is located, it has oscillat- 
ing character. When y « 2.5 A, fiy reaches its negative 
maximum. At this position, as schematized in Fig. [5Ji, 
the molecule experiences a negative electric field over al- 
most its whole size, whereas the average of F^ is almost 
zero and thus « 0. 

The knowledge gained above enables us now to search 
for the global minimum of the PES by changing the resid- 
ual degrees of freedom [l^ - translation of the molecule 
along the a;-direction and rotation around the z-axis. The 
QM/PPC approach makes it possible to scan a set of 
1500 different molecular configurations. The results are 
condensed in Fig. [3] The QM/PPC interaction energy 
(relative to the isolated molecule) is represented in Fig. 
[5^ as a function of the rotation angle 6 for the whole set 
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of x- and y-positions sampled over the ZnO(lOlO) surface 
unit cell. 
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FIG. 3; Adsorption scenario of 2P on the ZnO(lOlO) surface 
as computed by the QM/PPC method, a) Molecule-substrate 
interaction energy versus rotation angle 9 for all x- and y- 
positions counted. ^ = 0; long molecular axis along x- 
direction (see Fig la). In b) and c), the interaction energy 
is plotted versus the z- and j/-component, respectively, of the 
induced dipole moment 

The absolute minimum is found ai — 90° (long 2P 
axis lly). However, there is also a second minimum at 
9 = which corresponds to the one in Fig. [5^. The 
energy difference between the two minima is only 20 
meV and hence within the numerical error range. We 
conclude that the 2P molecule can be arranged on the 
ZnO(lOlO) surface in two different cross-aligned orienta- 
tions which makes the formation of a well ordered mono- 
layer rather questionable. 

The central question to be answered is about the mech- 
anism controlling the alignment of the molecule. The 
interplay between the surface electrostatic field, the in- 
duced dipole moment, and the interaction energy be- 
comes evident from Fig. [3)d and c. There is a distinct 
linear relation between the IS.E and /j,z. The energy is 
mimimized if and only if the dipole moment along z is 
maximized. On the other hand, as seen in Fig. |3j;, it 
holds /ij, « at the energy minimum, but there is no 
direct relation like in the case of fiz- 

The above findings become more transparent in an 
analytical model gathering the leading features of the 
molecule-substrate electrostatics. The energy of a 
molecule with zero static dipole but finite quadrupole 
moment Mg in a weak but non-uniform electric field 
(F, « 0) is O 
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where is the molecule's polarizability tensor. This 
expression is derived from perturbation theory: the first 
term represents the electrostatic interaction between the 
external non-uniform perturbing field and the unper- 
turbed molecule, the second one the induction energy ^IQ] 



accounting for the molecular polarization created by the 
field. Eq. ([ij is valid for a large class of planar, sym- 
metric oligomers characterized by vanishing off-diagonal 
elements of and aij. The ZnO(lOlO) surface peri- 
odic dipolar electric field seen by the molecule can be 
approximated by (see Fig 2c, d and Ref. [Tl|) 

Fy{y, z) w Ae-''"- cos(fcy) , F,{y, z) « -Ac-''"- sin(fcy) 

with k — 2TT/a so that its norm {F w Ae^''^) is indepen- 
dent on y and 
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dy dz 
Using that fii — auFi and inserting ([5]) in ([T]), we obtain 



AE « -B^z - Cfil 



a F^ 



with B = k{Myy — Mzz)/2azz and C — {ay 
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molecule 


Myy 






Ctzz 


B 


C 




a.u. 


a.u. 


a.u. 


a.u. 


eV/D 


eV/D^ 


2P 


-46.2 


-57.3 


217.4 


68.2 


0.558 


0.068 


6P 


-135.7 


-168.4 


1257.6 


190.5 


0.589 


0.062 


5A 


-82.6 


-101.4 


645.2 


115.2 


0.557 


0.084 


5PV 


-145.3 


-180.1 


2052.9 


207.8 


0.574 


0.090 



TABLE I: PBE quadrupole moments {Myy, M^z), polarizabil- 
ities {axxjCtzz) as well as B and C coefficients (see text) for 
difi'erent molecules. 



Table |T] compiles the values of the relevant param- 
eters for 2P and three other representative molecules 
(6P, 5A=pentacene, 5PV=penta-phenylene-vinylene) , all 
widely used in organic opto-electronics. Though the 
anisotropy of a can be quite significant, the term 
quadratic in /i^ in ([3]), originating from the induction 
energy, is negligible against the linear term (i.e., the 
quadrupolar contribution) for weak fields. Hence, the 
analytical model fully recovers the numerical results of 
Fig. a linear fit to the data in Fig. for 2P pro- 
vides a slope of -0.48 eV/D in very good agreement with 
the value expected from Table HI 

The analytical treatment suggests that the electro- 
static scenario found for 2P is of general validity. In- 
deed, the full numerical analysis confirms this for 6P, 6A, 
and 5PV. A global optimization for molecules of this size 
cannot be carried out using first-principles methods only 
and, usually, semi-empirical approximations have to be 
employed 17[. The QM/PPC approach is instead non- 
empirical because the only parameter (q) is fixed from 
a first-principles (PBE/PW) calculations. We thus were 
able to perform the same global scan of the PES as for 
2P. Fig. 0^ demonstrates that the linear relation between 



4 



interaction energy and vertical dipole moment is a com- 
mon feature for this class of non-polar molecules. Even, 
the slope of the curves is almost the same (-0.46 -0.58 
eV/D), consistent with the similar values of B in Tab. H] 
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FIG. 4: Interaction energy for 6P, 5A and 5PV on 
ZnO(lOlO) in QM/PPC. a) Linear relation between molecule- 
substrate interaction energy and induced vertical dipole mo- 
ment, b) Interaction energy versus rotation angle, only the 
minimum values from all x- and y-positions scanned are 
shown. Inset (color-online): Orientation of 6P at the global 
minimum shown on a colormap of the z-component of the 
surface electric field. 

Although the linear energy-dipole relation holds for all 
these molecules, the specific alignment on the substrate 
can be substantially different. This is documented in Fig. 

where the energy is plotted as a function of 9 at x- 
and y-positions with minimized energy. In contrast to 
2P, the PES of both 6P and 5PV exhibits deep global 
minima at = 0, clearly separated by 140 meV and 330 
meV, respectively, from other arrangements. Not only 
the orientation but also the lateral position of the ad- 
sorption site is uniquely defined, with y w 3.7 A and 
a; = for both 5PV and 6P. Thus, as illustrated in the 
inset for 6P, the energy is minimized when the long axis 
of the molecule matches with the lines of largest posi- 
tive Fz , where the electrostatic coupling and thus /i^ are 
maximized. The longer the molecule, the more stable 
the alignment. The PES of 5 A is instead less deep and 
structured. In contrast to 6P/5PV, no preferred orien- 
tation can be thus anticipated here. This finding can be 
rationalized by the fact that 5A has no carbon atoms ex- 
actly on the long molecular axis which can be most easily 
polarized by the electric field, as it is for 6P/5PV. 

In conclusion, we found that the periodic dipolar elec- 
tric field of the ZnO(lOlO) surface plays a key role in the 



adsorption of typical oligomers. When the molecules ex- 
hibit an axially oriented 7r-electron system, a well-defined 
molecular alignment, stabilized by energies larger than 
100 meV against reorientation, is established, as observed 
experimentally for 6P The electrostatic coupling is 
characterized by a linear relation between the molecule- 
substrate interaction energy and the induced vertical 
molecular dipole moment, which can be employed to pre- 
dict and/or to design the molecular orientation on the 
surface. Moreover, this dipole moment is directly associ- 
ated with workfunction changes and thus provides 
a tool for engineering the energy level adjustment of in- 
organic/organic hybrid structures Q- Finally, we note 
that the single-molecule adsorption described above will 
be perpetuated and will result in molecular assemblies 
reflecting the topology of the surface field. Although 
the induced dipole moment is modified by depolariza- 
tion effects [3, [l^ , this energy scale is certainly signif- 
icantly smaller than the electrostatic molecule-substrate 
coupling controlling the alignment on the surface. There- 
fore, we believe that our findings define a route towards 
the realization of highly-ordered self-assembled arrays of 
oligomer/polymers on ZnO(lOlO) and similar surfaces. 
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207441). 
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